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Abstract 

The phenomenon of Stochastic Resonance (SR) has been conclusively demonstrated in bistable 
potentials. However, SR in sinusoidal potentials have only recently been shown numerically to 
occur in terms of hysteresis loop area. We show that the occurrence of SR is not specific to 
sinusoidal potentials and can occur in periodic bistable potential, U{x) = |(cosx + cos 23;), as 
well. We further show that SR can occur even in a washboard potential, where hysteresis loops 
normally do not close because of average drift of particles. Upon correcting for the drift term, the 
closed hysteresis loop area (input energy loss) shows the usual SR peaking behaviour as also the 
signal-to-noise ratio in a limited domain in the high drive-frequency range. The occurrence of SR 
is attributed to the existence of effectively two dynamical states in the driven periodic sinusoidal 
and periodic bistable potentials. The same explanation holds also when the periodic potentials are 
tilted by a small constant slope. 
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I. INTRODUCTION 



The phenomenon of stochastic resonance (SR), initially put forth theoretically to ex 



plain the occurrence of ice ages, at a certain frequency, on earth 111, have been observed 

5^,0^. SR has been investigated 

a. 



experimentally in many physical and biological systems 

theoretically and experimentally with considerable interest over the past three decadesp, 
Its attraction lies in the seemingly counter-intuitive idea that by tuning noise level (exter- 
nally or internally) the response of a nonlinear system to a weak external periodic signal 
can be enhanced considerably; or a nonlinear system itself tunes the noise level in order to 
enhance a particular chosen signal. Moreover, the noise level at which the response peaks 
depends on the frequency of the input signal. As a consequence, it may have potential appli- 
cations in the detection of weak signals as well as in the selection of a signal of a particular 
frequency out of a host of signals of different frequencies js- 11 1. As a corollary, a biological 



system can tune noise level internally to select and enhance a desired signal 
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-ll^. Apart 



from potential practical applications, it offers considerable theoretical challenges. The oc- 



currence of SR in bistable systems 0, 

widelyQ, However, its occurrence in periodic potentials is still debatable [l6|, t 



15| has more or less been confirmed and re por ted fairly 

lough 



isj as 



there have been some investigations in monost able j 171] and periodic potential systems 
well. 

There have been some discussion on SR in periodic potentials in view of obvious prac- 
tical importance of these potentials, such as in describing motion of adatoms on crystal 



surfaces 



19|, superionic conductivity, RCSJ model of Josephson j unction |20N22|. etc. The 



conventional SR is considered to occur in a bistable system at a temperature when the sig- 
nal frequency matches the time rate of passage across the potential barrier of the nonlinear 



system 
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16|. Analogously, one would expect, for instance, the frequency dependent mobil- 



ity, as a response, to peak as a function of temperature, at a signal frequency corresponding 
to the mean passage rate across a potential barrier of the periodic potential. However, in- 
stead, the mobility shows monotonic behaviour around that frequency. It does show a peak 
with temperature only at a much higher frequency. This mobility peaking with temperature 
is termed as a kind of dynamical resonance, unlike the SR satisfying the conventional cri- 



terion of frequency matching, and merely confined to intrawell motion 
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-Il8|. This aspect 



has been examined more closely in a recent numerical work [23] on a sinusoidal potential 
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considering hysteresis loop area as response to the input signal F{t) = AFcosut. 

Hysteresis loop area (HLA), in the average position-forcing, (x—F), space, or equivalently, 
the input energy lost by the system to the environment per period of the external forcing 
F{t), is considered as an appropriate measure of SR. HLA has been considered as a quantifier 



of SR earlier too 



25[. However, recently HLA was found to show SR behaviour close to 



what was shown by the amplitude of the average particle-position variable x{t) in a bistable 



system 



26|-|29|. HLA not only includes information of the mean amplitude xq of x{t) but 



also its phase relationship with the external forcing or the input signal F(t). Its significance 



as an appropriate quantifier of SR has also been pointed out recently [301. 13l|. Moreover, in 
a periodic potential, it is very likely that the particle forays into wells far away, on either 
direction, from the initial well. It makes the position variable unsuitable for any meaningful 
quantification of SR. It is found that the HLA or input energy loss takes account of motion 
whether it is in a single well or spread over several wells. Also, HLA shows qualitatively 
similar behaviour as the frequency dependent mobility calculated using the linear response 
theory IGj. The HLA has recently been used as a quantifier of SR in a underdamped periodic 
potential system |23[]. 



In Ref. 



23| it was pointed out that the peaking of HLA as a function of teniperature 



at high frequencies ought not to be dismissed as merely a dynamical resonance 16| . The 
phenomenon can be seen as a result of transitions between two dynamical states of the 
damped driven particle trajectories in the highly nonlinear (sinusoidal) potential. The two 
and only two dynamical states exist and are distinguished by their trajectory amplitudes 
and their phase relationship with the field F{t). The relative stability of the two states 
changes as the temperature (noise strength) is varied. 

It is hard to prove analytically that an underdamped particle moving in a sinusoidal 
potential driven by a periodic force of frequency close to the natural frequency at the bottom 
of a well of the potential and subjected to fluctuating forces can have only two dynamical 
states of its trajectories as it was found numerically. Solution of only the simplest situation 
of the deterministic motion of a driven pendulum without friction. 



m 



— sinx) + Fq cos ut, 



given the initial condition (x(0), f (0)), can be obtained analytically and its phase trajectories 



drawn 



32l |. Notice that the periodic potentials Uo{x) = —sinx and Uo{x) = cosx are 



equivalent except for a phase difference of ±f . Also, the trajectories with initial conditions 
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(a;(0),0) and {2n — x{0),0) should be identical except that they are in opposite directions 
in case of closed orbits and in ±v sectors in case of travelling solutions. In other words, 
the two trajectories should have the same amplitude but different phase relationship with 
the external forcing. However, for different initial conditions (x(0),v(0)), they can have 
different trajectories and not necessarily confined to just two particular trajectories. But, 
when damping is present, it is easy to appreciate that if two trajectories have different phases 
they will necessarily have different amplitudes and hence different energy contents. This can 
be seen as follows. 




FIG. 1: The figure shows the periodic potential when F{t) = (curve O), F{t) = AF (curve ^4) 

and F{t) = -AF (curve B). 

Consider the two extreme cases of particle motion being (i) in-phasc, and (ii) completely 
out-of-phase with the external drive. Fig. 1 shows the periodic potential when F{t) — 
(curve O), F(t) = AF (curve A), and F(t) = -AF (curve B), where F(t) = AFcosa;^, so 
that the total potential V{x) — — sin a; — xAF cosut. Consider the extreme position PI of 
the particle on the curve A when F{t) — AF. In the next moment the F(t) will decrease 
and also the particle position moves to the left so that F{t) and x{t) are in phase. They 
being in phase implies that when F{t) = the particle is at the potential bottom P2 and 
when F{t) — —AF the particle is at P3, the extreme position on the left and is about to 
roll down the potential hill as F{t) begins to rise from — AF. Prom the figure, therefore. 
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one could notice that in this case of in-phase situation the particle always moves on the 
stiffer slope of the potential V{x). On the other hand, if we were to begin from the particle 
position P4 on A while F{t) = AF, the particle always moves lying on the gentler slope of 
V{x) in this completely out-of-phase case. Since the force experienced by the particle due to 
the potential in the two cases have different values one would naturally expect the damped 
particle to have different amplitudes of motion and hence have different energy losses due 
to frictional forces. The two dynamical states are thus distinct. 

At low temperatures, for given F{t), depending only on the initial conditions the system 



chooses one of the two dynamical states. These two states are quite stable [23j and no 
transition occurs between them. However, as the temperature is increased transition takes 
place between the dynamical states and, as a result, the relative population of these states 
changes. Thus, the mean amplitude of the trajectories and the overall phase when averaged 
over the initial conditions vary with temperature. The hysteresis loop area shows a peak at a 
temperature indicating stochastic resonance where the transition rate between the dynamical 
states acquires a particular value. It is also important to notice that well before the HLA 
peaks, the particle begins to surmount the potential barrier of V{x) and at resonance the 
motion is no longer confined to a single well of V{x)] the inter- well transitions become 
quite numerous. Of course, the inter-well transition rate is still quite low compared to the 



transition rate between the dynamical states 



23|. 



In the present work, we show that a driven underdamped particle exhibits stochastic 
resonance in (i) a periodic bistable potential, U{x) = |(cosx + cos 2a:), and also in (ii) 
washboard potentials (tilted sinusoidal as well as tilted periodic bistable potential). In this 
work we use HLA as a quantifier of SR but also supported by the signal-to- noise ratio (SNR). 

The appropriateness of the HLA as a quantifier of SR vis-a-vis the SNR has been discussed 
extensively in Ref . 30| . However, the SNR has been used in many earlier works to discuss SR. 
Without going into the relative merits of these quantifiers we present the SNR to complement 
extensive results that we show in terms of HLA. As has been pointed out 
do not peak at the same temperature. However, our main contention that SR is a distinct 
possibility in periodic potentials is reinforced by the results of SNR. 

In the case of bistable periodic potential too, effectively, two and only two dynamical 
states (one in-phase and one out-of-phase) of trajectories of a periodically driven under- 
damped particle are obtained. Since in a potential well of U{x) there exist two similar 



30| the quantifiers 
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FIG. 2: The periodic bistable potential U{x) = | (cos a; + cos 2x) is shown with two similar subwells 
in a well of the potential. 

subwells (Fig. 2) there are two (one in-phase and one out-of-phase) states in each subwell. 
The in-phase (out-of-phase) state in one subwell has the same amplitude and phase rela- 
tionship as the in-phase (out-of-phase) state in the other subwell. Therefore, energetically 
there exist only one in-phase and one out-of-phase states of trajectories. However, as time 
progresses, these trajectories can also be in wells and subwells other than the initial ones, the 
basins of attraction of the dynamical states can be identified in different wells and subwells. 
Consequently, the basins of attraction of the dynamical states become more complex than 
in case of sinusoidal potentials where only trajectories in different wells are distinguished. 

In the case of periodic bistable potentials HLA, as also the signal-to-noise ratio, shows a 
peaking behaviour, just as in the case of sinusoidal potential, as the temperature is varied. 
The non-monotonic behaviour of HLA as a function of temperature is related to the relative 
stability of the in-phase and out-of-phase states of trajectories. Moreover, the particle does 
make inter- well transitions also. 

The washboard potentials, that is a slanted sinusoidal potential and slanted bistable 
periodic potential, naturally yield finite particle drifts. As a consequence the hysteresis 
loops do not close. However, upon correcting for the drift factor, as explained in Sec. HI, 
the hysteresis loops close and the closed hysteresis loop areas could be made to agree with 
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the input energy loss. The HLA (or the input energy loss) so obtained again shows peaking 
behaviour as the temperature is varied. This is a clear signature of SR in the washboard 
potential. The occurrence of SR could again be explained in terms of transition between 
the two dynamical (one in-phasc and one out-of-phase) states of particle trajectories that 
are realized even in these washboard potentials. 

II. THE MODEL 

We consider motion of an underdamped particle along (i) a periodic bistable potential 
U{x) — |V^)(cosa; + cos2a:), (ii) a tilted sinusoidal potential Ui{x) — — T/osin/ca; — Fqx, and 
(iii) a tilted bistable potential U2{x) = U{x) — Fqx, where the tilt Fq represents a constant 
force. As mentioned earlier (Fig. 2) U (x) has two subwells of equal depth (bistable) in each 
periodic well of the potential. U {x) is symmetric about x = nn, where, n = =f1, =f2, ■ ■ ■ . The 
latter two potentials we shall refer to as the washboard potential and the bistable washboard 
potential, respectively, for -Po 7^ but small. In the bistable washboard potential U2{x) the 
two subwells exist but are no longer identical as they were in U (x) . 

A particle of mass m moving in a medium of friction coefficient 7 along a potential V{x) 
and driven by an external periodic forcing F{t) = AFcosut and subjected to a Gaussian 
white noise ^{t), is described here by the Langevin equation. 



The fluctuating forces C{t) satisfy: < ^{t) >= and < C(t)^(t') >= 25{t - t'). The 
temperature T is in units of the Boltzmann constant ks- We take V{x) = U{x), Ui{x), or 
U2{x) separately to discuss the nature of particle motion. 

The equation is written in dimensionless units by setting m — 1, Vq — 1, k — 1. The 
Langevin equation, with reduced variables denoted again now by the same symbols, corre- 
sponding to Eq. (2.1) is written as 



The noise variable, in the same symbol ^, satisfies exactly similar statistics as earlier. 




(2.1) 




(2.2) 
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III. NUMERICAL RESULTS 



We adopt the same numerical procedures as described in Ref. 23|]. The drive (signal) fre- 
quency (w = Y") is chosen to be close to but a little smaller than the natural frquency at the 
bottom of the wells of the potentials. However, they are not exactly equal to the respective 
natural frequencies (without damping). For the periodic bistable potential U{x), we take 
the period r equal to 4.8. The same optimum r is taken also for the bistable washboard 



potential U2{x). Similarly, we continue with r = 8 as in Ref. 23| for the washboard potential 
With the periods r of the external forcing F{t) = AF cos ut, we obtain the trajectories 



x{t) for given initial conditions (a;(0),f(0)) numerically |24l. l33| by solving the Langevin 
equation (2.2) and calculate the input energy, or work done by the field on the system, W, 
in a period r, as|34|: 

r-^HiMM,, (3,) 



a* 

where, the effective potential Ue{x{t),t) = V{x) — xF{t), and V{x) = U{x), Ui{x), or U2{x) 
as applicable. Therefore, 

rF{to+T) 

W{to,to + T) = - xdF = A, (3.2) 

where A is the magnitude of the HLA. The average input energy per period, W, averaged 
over an entire trajectory spanning A^i periods of F{t), is 

^ n=Ni 
^ n=0 

Typically, A^i, ranges between 10^ to 10^, as required. 

At very low temperatures W depends very strongly on the initial conditions (x(0), v{0)). 
This is because whether the trajectories are in the in-phase state (with a small phase differ- 
ence 4>i between x{t) and F{t)) or in the out-of-phase state (with a large phase difference 02) 
is determined by the initial condition. And, W depends on the kind of trajectory. Moreover, 
at such temperatures the particle remains trapped in the initial well of the potential and 
also no transition between in-phase and out-of-phase states of trajectories takes place in a 
well. However, as the temperature is gradually increased the intra-well transitions between 
the in-phase and out-of-phase states begin to take place and subsequently, inter-well tran- 
sitions also become frequent. Therefore, as the temperature is increased the dependence of 
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W on the initial conditions (x(0),v(0)) weakens. However, it is always sensible to ensemble 
average W over all possible initial conditions and obtain the average input energy per period 
<W>, which is also equal to the mean hysteresis loop area < A >, Eq. (3.3). 

As mentioned earlier, the response x{t) in the two dynamical phases of trajectories not 
only have different phase relationship with the forcing F{t) but the mean amplitude Xq of 
x{t) in the two cases are very different. Therefore, apart from the inherent stochastic nature 
of W ( Eq. 3.1) the values of W are different depending on whether, during the particular 
period (of F{t)) in question, the system is in the in-phase or in the out-of-phase state of 
trajectory or makes transition(s) between the states. A clear picture is revealed by the 
distribution P{W) of input energies at various temperatures. For sinusoidal potentials the 
evolution of P{W) with temperature furnishes important information on the occurrence of 



SR^ 



29|. Also, since the entire stretch of the trajectories consists of a mixture of in-phase 
(with phase and out-of-phase (with phase (j)2) states, the avera ge p hase lag must 
lie between (pi and 02. Naturally, is a function of temperature 35|, |36|]. The phase is 



obtained numerically 



231] by calculating the mean hysteresis loops < x{F{ti)) >, where, 

n=Ni 



x{F{U)) = —J2^iF{nr + U)), (3.4) 

n=0 

for all [0 < < r]. 

In the following subsections we present and discuss our numerical results separately for the 
three cases of underdamped particle motion in potentials U{x), Ui{x), and U2{x) driven by 
field F{t) at various temperatures T. We take the dimensionless friction coefficient 7 = 0.12 
and initial velocity f (0) = v{t = 0) = for all cases. The initial position x(0) = x{t = 0) 
are chosen at 99 equispaced points Xi,i = 1,2, ■■ ■ ,99 between the two consecutive peaks, 
e.g., [0 < a;j < 27r], for the potential U{x). Unless otherwise explicitely stated the amplitude 
AF of F{t) is taken equal to 0.2 and the tilt Fq = 0.1. 

In the case of washboard potential Ui{x) and the bistable washboard potential U2{x), 
with constant tilt Fq the particle does acquire a mean drift (of velocity v{Fq)) at elevated 
temperatures and consequently the hysteresis loops do not close. Therefore, a correction 
is required to make the loops close. In a period, F{t) changes from Fq — AF to Fq + AF 
and then back to Fq — AF in a cosinusoidal manner. During a period the particle moves 
on the average a distance of rU(Fo). Approximating the variation of F{t) to be linear the 
mean work done during a period r as a result of the mean particle displacement is equal 



to I X 2 X AF X Tv{Fo). Therefore, the required correction to the hysteresis loop area 
< A > due to the mean drift equals AFtv{Fo). This simple approximate correction closes 
the hysteresis loops thereby enabling us to calculate the HLA and equivalently the <W >. 




FIG. 3: Plot of particle trajectories x{t) (a) for U{x) at T = 0.0005, and (b) for Ui{x) at T = 0.001. 
F{t) (dashed line) is also included for comparison. The phase lags <^i ~ O.OStt and (1)2 ~ O.Stt (inset) 
for U{x) and <^i ~ O.OStt and ^2 ~ 0.67r (inset) for Ui{x). 

When the underdamped particle moves along the potentials U{x){— |(cosa; + cos2x)) 
and Ui{x){~ —sinx — Fqx) and driven by the external periodic field F{t) the trajectories 
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FIG. 4: The figure shows the basins of attraction for bistable potential at T = 0.0005. The 
dynamical states of the particle are represented by ip- (in-phase) and op- (out-of-phase) states 
with respect to applied periodic force. The wells and subwells in which the dynamical states exist 
are represented by the indices / and m within braces; e.g., op[0,l] indicates that the out-of-phase 
trajectory of the particle is in the (left) subwell-1 of the (initial) zeroth well of the potential. 

are essentially of just two kinds (states). The in-phase states correspond to trjectories 
x{t) which lag behind F{t) by a small phase (pi, whereas the out-of-phase states which lag 
behind F{t) by a large phase (p2- Of course, 0i and 02 are approximate average values 
which weakly depend on the potential and the temperature. At low temperatures (pi and (p2 
essentially maintain the same values at each period of F{t), Fig.3. The potential U{x) has 
two similar subwells (Fig. 2) in a well, the trajectories in either subwell are just the same two 
states. Here, these in-phase and out-of-phase states are identified by the symbols ip(/,m) 
and op(/,m), respectively. / indicates well number, for example, / = for the initial well, 
/ = — 1(+1) for the first well to the left(right) of the initial well, m indicates the subwell 
number: 1(2) for the left (right) subwell. These states are truly dynamical states. The basins 
of attraction of these states for the potential U (x) is given in Fig. 4, as an illustration at 
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temperature T — .0005. 
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FIG. 5: Plot of with x(0) for Ui{x) at T = 0.001 (a), and U{x) at T = 0.0005 (b). 

At temperature T = .0005 the two states are quite stable and no transition between 
them could be observed. Though the potential barrier between any two consecutive wells 
have almost the same value for U{x) and Ui{x), their well bottoms are quite dissimilar and 
hence the input energy per period W are very different, Fig. 5. For Ui{x), W are about 
0.1 (in-phase) and 1.15 (out-of- phase) , whereas for U{x) they are about 0.04 (in-phase) and 
0.32 (out-of-phase) . Naturally, transitions between the two states for U{x) occur at lower 
temperature than in case of Ui{x), Fig. 6. 

In the case of U{x) by the temperature T = 0.001 the indications of out-of-phasc state 
going to the in-phasc state, in a subwell, could be found and by T = 0.0015 a substantial 
fraction of the out-of-phase states have jumped to the in-phase state. At T = 0.002 all 
the out-of-phase states have made way to the in-phase state due to thermal fluctuations. 



(a) 
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FIG. 6: Plot of W with a;(0) for U{x) at T = 0.001 (a), Ui{x) at T = 0.004 (b). 

Therefore, close to T = 0.002 the system have the lowest average input energy < W > {or 
< A >) per period of F{t). And, from upward of T = 0.003, the in- phase states begin to 
jump to the out-of-phase state. At T — 0.005 even inter-subwell transitions could also be 
observed. The corresponding temperatures for the case of Ui{x) are much higher. 

For the system with potential Ui{x) the out-of-phase states begin going over to the in- 
phase state at around T = 0.004 and the process completes at about T = 0.006. The 
transition from the all-in-phase state to the out-of-phase state begins at about T = 0.009. 
Interestingly, however, the inter- well transitions for Ui{x) takes place at a much lower tem- 
perature (about T — 0.01) than in case of the system with potential U{x) which is only at 
about T = 0.23, Fig. 7. 

As stated earlier the periodic bistable potential U (x) have essentially two states because 
the two subwells are energetically identical. In the case of the bistable washboard potential 
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FIG. 7: Plot of x{t) shows an inter- well transition for Ui{x) at T = 0.01 (a). A similar transition 
is seen for U{x) at T = 0.23 (b). 

U2{x){— U{x) — Fox) the two subwells become dissimilar and the right subwell is energetically 
lower than the left one and we get two states corresponding to each subwell. However, the 
in-phase state in the left subwell becomes unstable for amphtude AF > 0.19 of the drive 
F{t). Therefore, in our case, we have only one in-phase state to consider in the left subwell. 

In case of U2{x), W is about 0.27 for the out-of-phasc state in the left subwell (henceforth 
called subwell-1) and W for the two states in the other subwell (subwell-2) are about 0.03 
(in-phase) and 0.37 (out-of-phase) at T = 0.0005. By the temperature T — 0.001 most of 
the (out-of-phase) states in subwell-1 go over to the states in subwell-2 and by T = 0.0015 
no states in subwell-1 survives. As the temperature is gradually increased the out-of-phase 
states in subwell-2 begin to jump over to the in-phase state in the same subwell and by 
T = 0.003 we only have in-phase states in subwell-2 and < W > acquires a minimum value. 
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By T — 0.009 the particles begin to leave the in-phase state for the out-of-phase state in 
subwell-2 and we get a mixture of the two states in the same subwell. As the temperature 
is increased further, the transitions back and forth to subwell-1 begin at about T — 0.015. 
Inter-well transitions begin, at a very slow rate, at about T = 0.08, which is much higher 
than the corresponding temperature for the washboard potential Ui{x) but lower than that 
for U{x). The presence of two subwells in U{x) and U2{x) makes the effective bottom of 
the wells flatter (smaller curvature) than in case of Ui{x) and hence the inter- well transition 
rates smaller (see Eq. (3.5) below). 

At T > 0.003 (T > 0.01) transitions from in-phase states to the out-of-phase states 
increase rapidly for U{x) {Ui{x)), and consequently so docs the <W >, leading ultimately 
to the SR condition. Similar is the situation for U2{x). 

B. Hysteresis loss and stochastic resonance 
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FIG. 8: Plot of {W) and SNR as a function of T for the washboard potential Ui{x). 

Figure 8 shows the variation of average input energy < W > or the average hysteresis 
loop area < A > and the signal-to-noise ratio (SNR) as a function of temperature for the 
washboard potential Ui{x). The same functions are plotted for the potentials U{x) and 
U2{x) in Fig. 9. <W > peaks at the temperatures T = 0.2, 0.04 and 0.08 for the potentials 
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FIG. 9: Plot of {W) and SNR as a function of T for U{x) (a) and U2{x) (b). 

Ui{x), U{x) and U2{x), with peak values of about 0.36, 0.16, and 0.146, respectively. These 
maximum values of < 1^ > lie between the W values of the in-phase and out-of-phase states 
for the corresponding potentials. There is, however, a remarkable difference between the 
nature of motion of particles at the temperature (Tsr) of maximum < W > for Ui{x) on 
one hand and U{x) and U2{x), on the other. 

As the temperature T goes through Tsr ~ 0.2 corresponding to maximum < W > , 
the particle feels the periodic nature of the potential Ui{x) as the inter- well transitions are 
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quite frequent, though not as frequent as the intra- well transitions between the in-phase 
to out-of-phase states. The particle motion is truly in a periodic potential implying the 
presence of stochastic resonance in the periodic potential Ui{x). On the other hand, the 
particle has not yet begun the inter-well transitions (in U{x)) or have just started (in U2{x)) 
as the temperature goes across the < W > maximum. Hence the particle does not feel the 
periodicity of the potentials but the presence of the two subwells. Therefore, the resonance 
effectively occurs in a single well with two subwells. However, the bistability in a well of 
the potentials is required for that to happen. Moreover, the existence of dynamical states 
of trajectories is necessary for the occurrence of < > maximum at that low temperature, 
TsR = 0.04 for U{x) and Tsr = 0.08 for U2{x). 

Consider the simple Smoluchowski limit of Kramers rate [s^] : 

fc = ^^exp(^), (3.5) 

where Uq and u^, respectively, are the curvatures at the bottom of the wells (subwells) 
and at the top of the barrier across the wells (subwells) of the potential Ui{x) {U{x) and 
U2{x)). This rate calculation shows that k^^ ~ 3450 across the potential barrier between 
two consecutive wells of Ui{x) (r=0.2, 7 = 0.12, = 1.686) and ~ 47 x 10^ and 
760 across the potential barrier between the two subwells in a well of the potentials U{x) 
(T=0.04, 7 = 0.12, Eb = 0.75) and U2{x) (T=0.08, 7 = 0.12, Eb = 0.618), respectively. 

Note that the periods r = 8.0,4.8 and 4.8 of the drive field F{t) taken, respectively, 
for the potentials Ui{x),U{x) and U2{x) were too small compared to the calculated k~^. 
Therefore, the observed SR cannot be considered as the conventional stochastic resonance 
following Ref.[l6|. However, the resonance is brought about by the existence of and transition 
between the two dynamical (in-phase and out-of-phase) states of trajectories. The rates of 
these transitions are of the order of the period r of the drive field. Viewed in this perspective 
of transition between the two dynamical states they, indeed, indicate SR. Moreover, the 
potential Ui{x), at T = Tsr a substantial number of inter-well transitions is observed. But 
in case of U{x) and U2{x), inter-well transitions are seldom observed at T = Tsr. However, 
inter-subwell transitions are comparable in number to the transitions between the dynamical 
states at Tsr. SR nature of these transitions is further supported by the behaviour of input 
energy distributions P(W) across the temperature (Tsr) of maximum < W >. 
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C. Input energy distribution and SR 



The input energy distributions P{W) have earher been used in the discussion of SR|23l. 
27H29|. P{W), in the present case are shown in Fig. 10. At low temperatures, say T = 
0.0005, P{W) shows the usual bimodal distribution for the potentials U{x) and Ui{x). The 
peaks occur exactly aX W = W values shown in Fig. 5, corresponding to the W in the in- 
phase and out-of-phase states. Fig. 10 also shows P(W) for the potential U2{x) exhibiting 
three peaks for the amplitude AF = 0.2 of F{t)- This is because, as mentioned earlier, the 
two subwells are not equivalent because of the finite average tilt Fq = 0.1 in the potential 
U2{x). On general considerations one would expect four peaks. The four peak P{W) appears 
only with AF < .19. The three peak P{W) instead occurs because of the instablity of the 
in-phase state in the subwell-1 for AF > .19. Here again the peaks are centered ai W = W 
values corresponding to the three dynamical states. The correspondence of the peaks of 
P{W) and the dynamical states is thus unambiguous at low temperatures. 

In Fig. 11 are shown P{W) at various elevated temperatures. As the temperature 
is gradually increased two important common features could be readily observed: (i) the 
variation of the strength of the P{W) peaks, and (ii) the intrusion of the in-phase peak into 
the negative W domain and the appearance of a long negative W tail of P{W). The analysis 
of these two features provides a better understanding of SR in the three potential systems. 

For the potentials Ui{x) and U{x) the higher energy peak (corresponding to the out-of- 
phase state) first diminishes, disappears and then reappears, as the temperature is gradually 
increased. For U2{x), the peak of P{W) corresponding to the out-of-phase state in the 
left subwell, i.e. op(0,l), first diasappears, and then the peak corresponding to op(0,2) 
too diasappears. As the temperature is increased further a broad peak (almost a plateau) 
appears roughly spanning the earlier two out-of-phase peaks. As the temperature is increased 
further the newly formed peaks begin merging with the sole in-phase peak (for AF = 0.2) 
and at the temperature Tsr the out-of-phase peak is left only as a receding shoulder leaving 
no distinguishable trace of either a hump or a plateau. However, in the process the strength 
of the in-phase P{W) peak also diminishes. This can be seen in either of the two ways: (i) 
by fitting the in-phase P{W) peak by a Gaussian and calculating its area and (ii) by finding 
the ratio of the number of points in the stroboscopic (Poincare) plots falling in the in-phase 
region to the total number of points. 
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FIG. 10: Plot of P{W) at different values of T for Ui{x) (a), U{x) (b), and U2{x) (c). For small 
T, the distribution is bimodal for Ui{x) and U{x). However, P{W) shows three peaks for U2{x). 
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FIG. 11: Plot of P{W) at different values of T for Ui{x) (a), U{x) (b), and U2{x) (c). At Tsr, 
P{W) is on the verge of losing its multimodal character; at higher T the distribution has a single 
peak structure. 
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The fraction (contribution) of in-phase states in the trajectory reduces from 1 (at T = 
Tmin ~ 0.02,0.003 and 0.003 for Ui{x), U{x) and U2{x), respectively, corresponding to 
the respective minimum of < >) gradually and becomes almost equal to 0.5 at the 
temperature Tsr of maximum <W > {Tsr ~ 0.2,0.04 and 0.08 for Ui{x), U{x) and U2{x), 
respectively). Fig. 12. At this temperature, it becomes hard to clearly distinguish the 
regions of in-phase and out-of-phase states of trajectories, just as in the hquid-gas system at 



the critical temperature. T 
of kinetic phase transitions 



le temperature of maximum < W > is said to fall in the region 



The second important feature of P{W) is its intrusion into the negative W region. As 
mentioned earlier it has two components: (i) a systematic broadening of the in-phase peak 
of P{W) and spilling over toW<0 region and (ii) the emergence of a long W < tail. The 
former happens at a temperature much lower than the temperature at which the W < 
tail begins to emerge. The phase lag (pi in case of in-phase is small {—4>i ~ O.Itt, O.OStt and 
O.OStt, respectively for the Ui{x), U{x) and U2{x) potentials). Also the in-phase P{W) peak 
is centered close to W > 0. As the temperature is increased the fluctuations in the value of 
(pi occur leading sometime to make (pi > 0. In other words, the thermal effect, makes x{t) 
occasionally lead the forcing F{t). Thus, a few individual HLAs acquire a sign opposite to 
what is observed in the usual case when causality is respected. Therefore, occasionally, W 
becomes negative and the in-phase peak of P{W) broadens into the W < region. This 
process continues and becomes more evident as the temperature is increased. 

Analyses of the trajectories x(t) and the accompanying hysteresis loops x{F{t)) reveal 
that the violation of causality (leading to W < 0) occurs during the transitions between 
the dynamical (in-phase and out-of-phase) states. In these situations the magnitude of 
negative W (hysteresis loop area with opposite sign) often becomes very large. The origin of 
negative tail of P{W) lies in these transitions between the dynamical phases. Since, before 
the temperature at which < W > becomes minimum at most one (out-of-phase to in-phase) 
transition occurs in the entire history of any trajectory spanning about 10^ periods of F{t) 
the negative tail of P{W) does not show up till the temperature Tmin of minimum < W > 
(Tmin ~ 0.02,0.003 and 0.003 for Ui{x), U{x) and U2{x), respectively). 

As the temperature is gradually increased further, the transitions, initially from the (all) 
in-phase to the out-of-phase states take place and the high energy peak (plateau) reappears. 
Thereafter, transitions in both directions become more and more frequent ultimately causing 
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FIG. 12: The figure shows the plot of the fraction (ratio) of time spent by the particle in its in- 

phasc state of trajectories as a function of temperature T, for the potentials U{x) {Tsr = 0.04), 
U2{x) {Tsr = 0.08), and i7i(x) {TgR = 0.2). 

the newly emerged high energy (out-of-phase) peak (plateau) to merge with the in-phase 
peak at the resonance temperature Tsr. The increase of negative tail of P{W) goes hand in 
hand with increasing T. At temperatures T > Tsr the transition regions dominate over the 
regions of in-phase and out-of-phase dynamical states. At these temperatures the increasing 
long negative tail and the ever broadening in-phase peak of P{W) bring down < W > from 
its maximum value at (Tsr). Note that on the average causality is always respected and 
< W > never becomes negative. 

D. Average amplitude and phase of hysteresis loops 

The average hysteresis loops < x{F{ti)) > are calculated using Eq. (3.4). Since the 
amplitude AF of the external drive F(t) — AFcos(^t) (=0.2) is much smaller than the 
periodic barrier heights (~ 2.0) between two consecutive wells, the amplitude xq of the 
average response < x(F(tj)) > is small and the relation < x{F{ti)) >= xocos(^t -|- 0), 
with t — riT + ti, is found to follow quite well for all three potentials at low temperatures. 
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FIG. 13: Plot of xo and -cj) (insets) with T for Ui{x) (a), for U{x) (b), and for U2{x)) (c). 

However, at higher temperatures the relation serves reasonably well for the potential Ui{x), 
and only approximately for U{x) and U2{x). 

The amplitude Xq and phase acquire their respective minima at the temperature T^j„ 
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of minimum < W > corresponding to the sole dynamical (in-phase) state, Fig. 13. At 
this temperature ~ (;/>!. As the temperatrure is gradually increased from T = Tmin xq 
peaks at T ~ Tsr. However, (p shows a monotonic behaviour. At these temperatures 
1^.1 < 151 < 1^.1 because the t.aiectone= co„.st of a n„xtu.e of in-„ha.e and out-ot-pha=e 
states. The variation of (/) is similar to what is reported earlier |35| and does not show a 
peak {36 1 . 



IV. DISCUSSION AND CONCLUSION 



The existence of two dynamical states in a sinusoidally driven (under) damped system is 
quite special to periodic potentials. The (quadratic) harmonic potentials or the (quartic) 
Landau potentials show only one (in-phase) state. The periodic nature of the potentials 
allows the particle to explore (in space) regions of high nonlinearity. The high nonlinearity of 
periodic potentials appears to be responsible for the occurrence of these two states (especially 
the additional high amplitude out-of-phase state) in a well in a sinusoidal potential or in a 
subwell in the bistable periodic or washboard potential. 

The investigation of dependence of the existence of two dynamical states on the amplitude 
of the drive field could also be of interest in periodic potentials. As observed earlier, the 
in-phase state in the left subwell of the potential U2{x) disappears for AF > 0.19. However, 
if AF is chosen to be too small the particle may not have the opportunity to explore the 
highly nonlinear regions of the potential and the out-of-phase trajectories may disappear. 
Therefore, the choice of AF too is crucial for the study of SR in these potentials. 

In the conventional SR it is the bistability of the potential that plays a crucial role. In 
the present periodic (or washboard) potential case SR is brought about by the bistability 
(or multistability) of the states of trajectories in a well (subwell) together with the adjacent 
wells (subwells) of the potential. Moreover, if the system is initially prepared in a given 
well of a double-well potential the conventional SR occurs when the probability distribution 
of particles becomes equal in both the wells in a finite (small) time (typically, of the order 
of the period of the drive field). In other words, the rate of passages across the potential 
barrier between the two wells becomes equal. Exactly similar is the condition in the periodic 
potential cases considered in the present work. The ratio of time spent by the particle in 
either of the two states of trajectories to the total time reaches 0.5 as the temperature rises 
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to TsR, Fig. 12. 

In order to observe SR in the 'periodic' potentials discussed above it is necessary to choose 
the period r of the drive field judiciously. SR occurs in a narrow 'window' of frequency u 
{u = ^): approximately, [7.8 < r < 9.8] for the potential Ui{x) and [4.8 < r < 10.6] for 
U{x) and [4.0 < r < 9.8] for the potential U2{x). The lower limit of the r-window is sharp. 
However, the upper limit is not well demarcated. The Tsr increases monotonocally with r. 
Here the upper limit of the r-window is arbitrarily put so that Tsr — 0.7. Tsr — 0.7 is a large 
temperature and is comparable to the potential barrier ~ 2 between two adjacent wells of the 
potentials. Consequently, < W > (T) peak becomes broader with increasing r. However, 
the numerical results for the potentials U{x) and U2{x) are expected to improve if r is 
chosen larger than the value 4.8. Indeed, the preliminary results show that when r is chosen 
somewhat larger than 4.8 even frequent inter-well transitions could be observed around Tsr 
for the potentials U{x) and U2{x). This effectively addresses the earlier limitation of particle 
motion being confined to a single well. This provides further support to the thesis of the 
presence of SR in these 'periodic' potentials. 
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